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The evolution of the centre-of-mass wave-function for a mesoscopic particle according to the 
Schrodinger-Newton equation can be approximated by a harmonic potential, if the wave-function is 
narrow compared to the size of the mesoscopic particle. It was noticed by Cohn et al. [Phys. Rev. 

A, 93, 062102 (2016)] that, in the regime where self-gravitational effects are weak, intermediate 
and wider wave-functions may be approximated by a harmonic potential as well, but with a width 
dependent coupling, leading to a time evolution that is determined only by a differential equation 
for the width of a Gaussian wave-function as a single parameter. Such an approximation results 
in considerably less computational effort in order to predict the self-gravitational effects on the 
wave-function dynamics. Here, we provide an alternative approach to this kind of approximation, 
including a rigorous derivation of the equations of motion for an initially Gaussian wave packet, 
under the assumption that its shape is conserved. Our result deviates to some degree from the 
result by Cohn et al., specifically in the limit, of wide wave-functions. 


I. INTRODUCTION 


the one-particle 


Schrodinger-Newton equation II, 


As long as there is no agreement on what the correct 
theory is that unifies quantum physics and gravitation, 
there is no unambiguous prediction for the gravitational 
interaction of quantum matter. While most physicists 
prefer the idea that the gravitational field itself must ex¬ 
hibit quantum properties at the microscopic regime, a 
fundamentally semi-classical coupling of quantum mat¬ 
ter to classical space-time has been proposed as an al¬ 
ternative M- A natural way to source the curvature 
of classical space-time by quantum matter is provided by 
the senri-classical Einstein equations 0,1, 

R ltv -\g IM ,R=^{*\f IM ,\V), ( 1 ) 

where the stress-energy-tensor of the classical Einstein 
equations is replaced by the expectation value of the 
corresponding quantum operator. While some contro¬ 
versies exist about the relevance of this approach 0,1, 
the arguments against it are inconclusive 0-(a|. The ul¬ 
timate decision whether or not the gravitational field 
must be quantised, therefore, seems to be up to experi¬ 
ments 15. 3, Q- 

Interestingly, the mass regime where such experiments 
become feasible is not that far from what is currently 
possible in the laboratory. This is because eq. © yields a 
nonlinear Newtonian gravitational self-interaction in the 
nonrelativistic Schrodinger equation. Such a nonlinearity 
has observable consequences, despite the weakness of the 
gravitational interaction. 

These self-gravitational effects were first modelled by 
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U(f,i-). (2) 


The self-gravitational interaction represented by the non¬ 
linear term predicts inhibitions of the free spreading of a 
wave-function for masses of about 10 10 u and beyond 
Experiments aiming at testing this behaviour could be 
realistic, at least in space, where a free spreading over 
longer time scales is feasible. 

Evidently, the massive systems needed for such an ex¬ 
periment cannot be elementary particles. Rather will 
they consist of a large number of constituents. One 
may, therefore, ask in what sense eq. © describes such 
a many-body system. The many-body dynamics follow¬ 
ing from the hypothesis of a fundamentally senri-classical 
gravity © will be discussed in some detail in section El 
One finds that—quite contrary to the situation where it 
is considered as a Hartree approximation 0 —eq. © 
approximately describes the centre of mass of a complex 
system, provided that its centre of mass is sufficiently de¬ 
localised in order to be allowed to disregard its internal 
structure. 

On the other hand, if the centre-of-mass wave-function 
is narrow compared to the size of the system, the result¬ 
ing self-gravitational potential is a function of the first 
and second moments, (r) and ( r 2 ), and quadratic in r. 
A potential of this form will leave the shape of a Gaus¬ 
sian wave-function intact during its evolution. Contrary 
to eq. ©, where for a numerical analysis the full wave- 
function must be simulated, the centre-of-mass dynam¬ 
ics can then be described by the evolution of the wave- 
function width as a single parameter. 

Colin et. al jTlj make use of this, proposing to approx¬ 
imate the Schrodinger-Newton equation by a quadratic 





2 


potential also in the case where the full dynamics are 
given by eq. ©• Their approximation is valid for weak 
gravitational potentials and substantially simplifies nu¬ 
merical calculations. It has been employed in order to 
provide a crucial contribution towards an experimental 
test by comparing the self-gravitational inhibitions of dis¬ 
persion to other decoherence sources. 

The wave-function dependent quadratic potential they 
use was, however, inferred from a comparison of ener¬ 
gies, rather than derived directly from the Schrodinger- 
Newton equation. Here, after reviewing the many-body 
and centre of mass equations in section im we will pro¬ 
vide such a direct derivation in section uni The result 
will be compared with the approximation used by Colin 
et al. [Ill, [jjj in section El where we find some devia¬ 
tions. We will conclude with a brief discussion. 


II. THE MANY-BODY 
SCHRODINGER-NEWTON EQUATION 


In a condensate of many constituents which collectively 
occupy the same state, eq. © appears as a Hartree ap¬ 
proximation and can be derived in a similar fashion as the 
Gross-Pitaevskii equation fl3l - [l5j |; it is then, for exam¬ 
ple, a model for the self-gravitation of Bosonic stars M- 
In this context, its origin is the mutual Newtonian in¬ 
teraction of the constituents and it provides an effective 
description of the factors of a tensor product state of the 
condensate, regardless whether eq. © is fundamentally 
correct or not. This is not the point of view taken here. 

On the other hand, if one starts with eq. © as a hy¬ 
pothesis about the gravitational interaction of quantum 
matter, applies it to a Klein-Gordon or Dirac field, takes 
the nonrelativistic limit, and restricts oneself to the one- 
particle sector then eq. © follows immediately as the 
fundamental dynamical equation for a single , nonrela¬ 
tivistic quantum particle 17, ll8j■ Since in the nonrela¬ 
tivistic limit there is no interaction between the sectors 
of different particle number, one can choose a restric¬ 
tion to the V-particle sector instead, for arbitrary N . In 
contrast to the Hartree approximation which is already 
an effective description for N constituents, one then ob¬ 
tains a different equation for the dynamics of the full 
IV-particle wave-function EM!: 

.. d . ( A n 2 _ 2 

v i —1 
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r i and m,i are the coordinates and masses of the N con¬ 
stituents. The gravitational potential V g contains both 
the self-gravitational interaction of each constituent and 
the mutual gravitational interactions between different 
constituents. 

Due to the nonlinear gravitational term, eq. (I3al) does 
not strictly separate into centre of mass and relative co¬ 
ordinates. Nonetheless, one can apply an appropriate, 
Born-Oppenheimer-type, approximation scheme, making 
use of the fact that the gravitational interaction is much 
weaker than intramolecular forces, which are making up 
the internal structure of the system modelled by a mass 
density distribution p c . Precisely speaking, if the wave- 
function separates into centre of mass and relative coor¬ 
dinates, ^(t; ri, • • • ,rjv) = if(t,r)\(r u ,riv-i), then 
p c is given by the relative wave-function 20]: 


N—l 

Pci r) =J2 mi 
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The centre-of-mass wave-function then satisfies 0] 
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Here r is the centre-of-mass coordinate, if the centre-of- 
mass wave-function, and p c the mass density of the many- 
body system with respect to its centre of mass. The 
function I Pc is proportional to the gravitational energy 
between a mass distribution p c and a copy of the same 
mass distribution shifted by a distance d. Equation © is 
the Schrodinger-Newton equation for the centre of mass 
of a complex system of many constituents. This is the 
equation we are interested in. 

If the system is modelled by a homogeneous, spherical 
mass distribution, the potential I Pc takes the form [2lj | 


dpf ere (d) = ^ 

x {!- 2(^) 2 + |(^) 3 - i (^) 5 ( d <2 R ), (6) 

(d>2R). 

One can see immediately that in situations where the 
mass distribution is point-like compared to the wave- 
function, eq. © is well-approximated by eq. ©■ 

In the opposite case, where the spatial wave-function 
can be considered narrow, and the mass is distributed 
homogeneously over a sphere of radius R , a good ap¬ 
proximation to V g is given by [2(1 [22[ 


V ™= const. 

+ ^^( r 2 - 2r-(tp\ r\if) + (if\r 2 \ ^}) - (7) 
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For a mesoscopic, homogeneous sphere one finds w| N = 
Gm/R 3 . The nonlinearity is still present, encoded in the 
expectation values. Based on this approximation, tests 
of the Schrodinger-Newton equation have been proposed 
also with optomechanical experiments [Hl - l24l |. 

Since the potential 0 depends only on the first and 
second moments and is quadratic in r, it can be shown 
to leave the_shape of an initially Gaussian wave packet 
unchanged [H], [H, HH] . Therefore, the full (as far as the 
narrow wave-function approximation is valid) dynamics 
for a Gaussian wave packet are given by the time evolu¬ 
tion of the width \J ( r 2 ) of the wave packet. 


III. EVOLUTION OF A SELF-GRAVITATING 
GAUSSIAN WAVE-FUNCTION 

For wider wave-functions, the gravitational potential 
deviates from the quadratic form 0, and the wave- 
function shape will change over time. If the self- 
gravitational interaction is weak (i. e. the kinetic energy 
of the mesoscopic particle is much larger than the gravi¬ 
tational energy), one may, however, assume as a first or¬ 
der approximation that the shape of the wave-function 
stays a Gaussian, and self-gravitation only affects its 
width Pill. 

In the case of a spherically symmetric wave-function, 
the expectation values of radial momentum and distance 
vanish exactly, (r) = (p) = 0. The equations of motion 
for the second moment u(t) = ( r 2 )(t ) follow straightfor¬ 
wardly from the Schrodinger-Newton equation ([5)1 . One 
finds 

J^X*) = ^;(r-p + p-r)(t) (8a) 

|(p 2 )(t) = -( p .(W g ) + (W g )-p) 

= - m ^< V ®)(*) ( 8b ) 

J^(r-p + p-r)(i) = ^-(p 2 )(t) - <r-Wg)(i), (8c) 

where the second equality in eq. (I8bl) is shown in ref. [23[ . 
If now one applies the aforementioned approximation, as¬ 
suming that the gravitational potential V g is simply a 
function of the width u(t) of a Gaussian state, this yields 

u (t) = -3wg N f(u(t)) u(t) (9a) 


For a homogeneous sphere mass distribution, the func¬ 
tion II takes the form 


n(z) = 


6 T 2_ 3-4 , 21 5_ 9. 7 ^<11 

t-us 2 ^ ' 20 ^ 70*^ — - L / ? 

\x (x > 1) . 


(9d) 


Note that the dependence on mass and the gravitational 
constant are completely encoded in wsn, and the function 
/ depends only on the wave-function width u, as well as 
the radius R of the mesoscopic particle. For the homo¬ 
geneous sphere, the integral can be solved analytically, 
resulting in 


/(») = erf (Jfj + 

( 7 324 - 162u - 35 u 4 + 70u 5 

V 2 70 u 4 

where erf is the Gauss error function and u is expressed 
in units of R 2 . One can easily check that this function 
has the appropriate limits 



lim f(u) = 1 and lim f(u) = 0. (11) 

u—t 0 u —^oo 


A. Initial conditions 

For an initial real valued Gaussian of width (r 2 )o = tto, 

' t(f=0 ’ , ' ) = (db) (-£;)■ (12) 

the differential equation (19al) must be solved with initial 
conditions 


u( 0) = u 0 

(13a) 

“o 

II 

o 

(13b) 

9h 2 

m(0) = 2?n 2 uq ^9(no)n 0 

(13c) 

with the function 
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For the homogeneous sphere, I\ is 



with (see the Appendix for a more detailed derivation) 


f(u) 
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f — -^-X^ — -^-X^ (x 1) 
K(x) = { 2 ^ 20 X 14 X 5 

(*>!)> 

which yields (again for units with R = 1) 
g(u) = erf 1 ' 


x | -u - 3 + 


486 + 105u 3 — 70u 4 _ 3/t 
105u 3 


(15) 


(16) 
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(a) f(u) (b) f (u) 




FIG. 1. Comparison of the function f(u) defined in eq. (I9al) (solid curve) and the function f(u) defined in eq. (1201) (dashed 
curve), u is in units of R 2 . (a) for small u. (b) for large u (logarithmic scale). 


The dynamics depend, in general, on the radius R of the 
mesoscopic sphere, its mass m, and the initial width uq 
of the centre-of-mass wave-function. However, the de¬ 
pendence on the latter is only a dependence on the ratio 
uq/R for the functions / and g , while wsn is a function 
of mass density only. For a given mass density, if units 
are chosen such that R = 1, the only true dependence on 
mass comes through the initial condition (II3d) . 

While for a precise numerical analysis of the equa¬ 
tions m and ([5]) the full wave-function must be simulated 
on a grid of many points, for this approximation it is suf¬ 
ficient to determine the evolution of u(t) = (r 2 ). This 
significantly simplifies numerical estimations of the mag¬ 
nitude of effects, and may even allow for some analytical 
results which cannot be seen from the full equation. 


IV. COMPARISON WITH PREVIOUS RESULTS 


Equation (KJal) also follows, if one starts with a har¬ 
monic oscillator ansatz with width dependent coupling, 

m^(t,r)= (-^-S/ 2 + ^-r^(t,r). (17) 

The authors of EMI employ an energy argument in 
order to determine an approximation for k(u). To follow 
this argument, first note that equation CB possesses the 
conserved energy 

E =^— /"d 3 x |V'0(f,r)| 2 + U(u(t)) , (18) 

2m J 


if the potential U satisfies 


dU k{u) 
drt 2 


(19) 


For the standard harmonic oscillator potential, where k 
is a constant, U is simply the potential energy of the 
oscillator. The potential U is then chosen such that 


U(u) = — G/pP here (- v /u), with J®P here defined in eq. ©, 
which leads to the functions 


f(u) 

d(u) 


1 -§Vu+ 


. 2v / u F 

1 - TeVu- 
. Vv? 


(u < 4) 
(u > 4) 

(u < 4) 
(u > 4) ’ 


(20a) 

(20b) 


rather than f{u) and g[u) in equations (liial) and (I13cl) . 

Intuitively speaking, the mesoscopic particle moves in 
a time-dependent harmonic potential, whose coupling 
constant is tuned in such a way that the potential en¬ 
ergy of the particle in the trap matches the gravitational 
potential energy between two particles, one sitting at the 
centre of the harmonic potential and the other one sitting 
at sj (r 2 ). 

From this intuitive picture it is already evident what 
could go wrong with such an approximation, if applied to 
the self-gravitational potential in the Sclirodinger-New- 
ton equation. For the case of a narrow wave-function, the 
gravitational potential within the mesoscopic particle is, 
in fact, a harmonic potential, as is also indicated by the 
approximation ©• However, in the wide wave-function 
case, where the mass density is well approximated by a 
delta distribution, the gravitational potential for a Gaus¬ 
sian wave-function with (r 2 ) = u is 



( 21 ) 


While this can be approximated by a quadratic function 
for small r, one obtains V g ~ 1/r for large r. There¬ 
fore, for the full Schrodinger-Newton equation the part 
of the wave-function at positions r 2 u will not see a 
quadratic potential, not even as an approximation, but 
a Coulomb-like potential. Hence, from the evolution ac¬ 
cording to the full Schrodinger-Newton equation one ex¬ 
pects a weaker response to gravity than the approxima¬ 
tion by a harmonic potential according to equations 113 
and m suggests. 
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(a) u(t) 



t (hours) 


(b) u(t) 



t (hours) 


FIG. 2. Solution of the equations of motion for a wide and intermediate wave-function. The dotted curve shows the free 
solution in absence of self-gravity, the solid curve is the solution using f(u) and g(u) as defined in eqs. (I9al) and (113cl) . and 
the dashed curve uses the approximation according to eq. ((20)). The mass is chosen to be m = 10 10 u for an Osmium sphere 
(p = 22.6 g/cm 3 ), corresponding to a diameter of approx. 100 nm. u is in units of R 2 . (a) wide wave-function, uq = (10f?) 2 . 
(b) intermediate wave-function, uo = R 2 ■ 


In Fig. Q] the functions f{u) and f(u) are plotted. 
There is a clear deviation in the intermediate regime, 
where u ss R 2 . For wide wave-functions, the deviation 
decreases in magnitude. Nevertheless, the functions also 
differ in their limiting behaviour for u —> oo. While 
f(u) approaches zero like ^3/(16^u 3 ), f(u) falls off like 
1 jyj 4 u 3 . Fig. [2] shows how this can lead to a qualitatively 
different behaviour in the equations of motion. The plot¬ 
ted solution corresponds to an Osmium sphere of about 
100 nm diameter. For a wide wave-function both solu¬ 
tions show at least a similar tendency to slightly decrease 
in width and oscillate, even though with a significantly 
lower frequency for the approximation with /(it), g(u) 
compared to f(u), g{u). In the intermediate regime, 
on the other hand, the solution according to eqs. (l9al) 
spreads only a little slower than the free solution, while 
the solution according to eqs. © and © is bound close 
to its initial width. 

One can use the equation (113cl) for the initial accel¬ 
eration, in order to define a critical mass (by the con¬ 
dition ii(0) = 0) which can serve as a rough estimate 
of whether the initial state will decrease or increase in 
width upon free evolution. Comparison of the functions 
g(u ) and g(u) then shows that, for wider wave-functions 
with u > (10i?) 2 , the approximation chosen in ref. flTIlT^ 
underestimates this mass value by an approximate 20%. 


V. DISCUSSION 

We have shown how, under the assumption of a weak 
self-gravitational effect that leaves the shape of an ini¬ 
tially Gaussian wave packet intact, the free quantum 
evolution according to the Schrodinger-Newton equation 
can be approximated by a single parameter differential 
equation for the square of the wave-function width, (r 2 ). 
The comparison with the previous approximation scheme 
used by Colin et al. [lj, E] for the very same purpose 


revealed deviations for wider wave-functions. While their 
results made use of an approximation of the gravitational 
potential by a quasi-harmonic interaction, which was de¬ 
termined by a rather ad hoc energy argument, here we 
provided a systematic derivation, based on the aforemen¬ 
tioned approximation. 

For the purpose of ref. Et where self-gravitating ef¬ 
fects were compared to decoherence sources, both ap¬ 
proximation schemes can serve as sufficiently good order 
of magnitude estimates. The behaviour in Fig. Or shows, 
however, that large deviations from the real evolution 
according to the full equation are possible. 

The alternative approximation scheme presented in 
this work is in no event inferior to the previously pro¬ 
posed one by Colin et al., as far as numerical efforts are 
concerned, since the functions f{u) and g(u) are known 
analytically. This scheme can also be easily adapted to 
mass density distributions other than the homogeneous 
sphere. 
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Appendix A: Derivation of (V g )(u) and (r ■ W 9 )(u) 

We want to derive (V g ) with V g as defined in eq. ( 0 , 
for the homogeneous sphere potential ([6]) . We write u = 
(r 2 ). The absolute value squared of the wave-function is 

w 5 <“> = (i) exp (y?r)' <A1) 
















First define a = 3R 2 /(2u) and 


with 
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§x 2 - ^x 4 + jqX 5 - ^a ; 7 (x < 1) , 
\x (X > 1) . 
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Then with appropriate substitutions we have 
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Further substituting 2 by x, and finally p = a + b, q = 
\a — b\, yields 


(V g )(a) 


512 Gm 2 a 3 

nR J 0 
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xe -4 a (a 2 +b 2 ) ah / dxxi(x) 
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(p 2 -q 2 )(J(p)^J(q)) (A5) 


In order to obtain (r • X7V g )(u) we have to repeat essen¬ 
tially the same calculation with the replacement 

i(x) —> ( a 2 — b 2 + x 2 ')——— , (A7) 

2x 

where the a 2 and b 2 terms cancel for symmetry reasons. 
This yields the same final relation 
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but with 
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(A9) 


[1] C. Mpller, in Colloques Internationaux CNRS, Vol. 91, 
edited by A. Lichnerowicz and M.-A. Tonnelat (CNRS, 
Paris, 1962). 

[2] L. Rosenfeld, Nucl. Phys. 40, 353 (1963) 

[3] S. Carlip, Class. Quantum Grav. 25, 154010 (2008) 
arXiv:0803.3456 [gr-qc] 

[4] K. Eppley and E. Hannah, Found. Phys. 7, 51 (1977) 

[5] D. N. Page and C. D. Geilker, 
Phys. Rev. Lett. 47, 979 (1981) 

[6] J. Mattingly, in Einstein Studies Volume 11. The Uni¬ 
verse of General Relativity , Einstein Studies, edited by 
A. J. Kox and J. Eisenstaedt (Birkhauser, Boston, 2005) 
Chap. 17, pp. 327-338. 


[71 J. Mattingly, Phys. Rev. D 73, 064025 (2006) 
arXiv:gr-qc/0601127 


[8] M. Albers, C. Kiefer, 

and 

M. R.egi- 

natto, Phys. Rev. D 78, 064051 (2008) 

arXiv:0802.1978 [gr-qc] 

[9] D. Giulini and 

A. 

Grofiardt, 

Class. Quantum Grav. 28, 195026 (2011) 


arXiv: 1105.1921 [gr-qc] 



[10] S. L. Adler, J. 

Phys. A 40, 755 (2007) 

arXiv:quant-ph/0610255 
[11] S. Colin, T. Durt, 

and 

R. Willox, 

Phys. Rev. A 93, 062102 (2016) 

i 


arXiv: 1402.5653 [quant-ph] 

[12] S. Colin, T. Durt, 

and 

R. Willox, 


Class. Quantum Grav. 31, 245003 (2014) 


arXiv: 1403.2982 [quant-ph] 

[13] E. P. Gross, II Nuovo Cimento 20, 454 (1961) 

[14] L. P. Pitaevskii, Sov. Phys. JETP 13, 451 (1961) 

[15] C. Bardos, L. Erdos, F. Golse, N. Mauser, and H.-T. 

Yau, C. R. Acad. Sci. Paris 334, 515 (2002) 

[16] R. Ruflini and S. Bonazzola, 

Phys. Rev. 187, 1767 (1969) 

[17] D. Giulini and A. Grofiardt, 

Class. Quantum Grav. 29, 215010 (2012) 

arXiv: 1206.4250 [gr-qc] 

[18] M. Bahrami, A. Grofiardt, S. Donadi, and 

A. Bassi, New J. Phys. 16, 115007 (2014) 

arXiv: 1407.4370 [quant-pli] 

[19] L. Diosi, Phys. Lett. A 105, 199 (1984) 

[20] D. Giulini and A. Grofiardt, 

New J. Phys. 16, 075005 (2014) 

arXiv: 1404.0624 [gr-qc] 

[21] H. I we, Z. Phys. A 304, 347 (1982) 

[22] H. Yang, H. Miao, D.-S. Lee, B. Helou, and 

Y. Chen, Phys. Rev. Lett. 110, 170401 (2013) 
arXiv: 1210.0457 [gr-qc] 

[23] A. Grofiardt, J. Bateman, H. Ulbricht, and A. Bassi, 

“Effects of Newtonian gravitational self-interaction in harmonically t 
(2015), arXiv:1510.01262 [quant-ph], 

arXiv: 1510.01262 [quant-ph], 

[24] A. Grofiardt, J. Bateman, H. Ulbricht, and 

A. Bassi, Phys. Rev. D 93, 096003 (2016) 

arXiv: 1510.01696 [quant-ph] 











